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1 Introduction 

The Traveling Salesman Problem (TSP) is the problem of finding a least-cost sequence in which 
to visit a set of cities, starting and ending at the same city, and in such a way that each city 
is visited exactly once. Interest in this problem cuts across several fields of study in Business, 
Economics, Engineering, and Computer Science, due partly to the wide applicability of the problem 
in practice (see [7] ) , and partly to the central place the problem occupies with respect to theoretical 
developments in the areas of Combinatorial Optimization (see [8]) and Computational Complexity 
(see [5] or [10] , among others). The seminal Operations Research work on the problem is that of 
(|H). A linear programming (LP) formulation was proposed in [3]. That formulation has 0(n 9 ) 
variables and (3(n 8 ) constraints (where n is the number of cities). 

In this paper we propose a new linear programming formulation of the Traveling Salesman 
Problem. The modeling approach used in the paper is similar to that of [3] . However, the proposed 
model is an order of size smaller, having 0(n 8 ) variables and 0(n 7 ) constraints. A small numerical 
experimentation we conducted shows that the computational times for the linear program proposed 
in this paper are several orders of magnitude smaller than those for the model in [3]. 

The plan of the paper is as follows. The proposed linear programming formulation is developed 
in section 2. The numerical experimentation is discussed in section 3. Conclusions are discussed in 
section 4. 



Notation 1 The following notation will be used throughout the rest of this paper: 



1. The number of cities is denoted n; 

2. The set of cities is denoted N := {1, 2, . . . , n}; 

3. The set of all TSP tours is denoted A; 

4. The cost of travel from city i to city j (j ^ i) is denoted 
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The set of real numbers is denoted by TZ; 
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For two column vectors a and b 




T ,b T ) T will be written as "(a, b)" (where 



(•) T denotes the transpose of (•)), except for where that causes ambiguity; 

7. The i th component of a column vector a is denoted a«; 

8. The notation "0" denotes a column vector of comfortable size that has every entry equal to 



9. The notation "1" denotes a column vector of comfortable size that has every entry equal to 



10. The convex hull of (•) is denoted Conv(-); 

11. The set of extreme points of (•) is denoted Ext(-). 

2 Development of the Formulation 

Assumption 2 We assume without loss of generality (w.l.o.g.) that: 

1. The number of cities is greater than 5 (i.e., n > 6); 

2. City 1 is the starting and ending point of travel; 

3. All vectors of variables are column vectors. 

Notation 3 



2. The set of "times-of-travel" to the cities in M is denoted T := {l,...,n — 1}; (i.e., if the 
time-of-travel to city i G M is r (r G T), then i is the r th city to be visited after city 1); 

3. The 0/1 binary that indicates whether city i is visited at time(- of -travel) r is denoted Wi r ; 
(wir = 1 iff city i is visited at time r); 

Theorem 4 There exists a one-to-one correspondence between the TSP tours and the perfect 
matching s of M and T; That is, there exists a one-to-one correspondence between TSP tours and 
points of 



Proof. Trivial. ■ 

The basic idea of our approach in this paper is to model the TSP as an optimization problem 
over a reformulation of the linear programming relaxation, 



0: 
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M:= A\{l} = {2,...,n} 




(1) 
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of Wi. Hence, we do not deal with the TSP Polytope per se (see [61 pp. 256-261]). Hence, results 
pertaining to descriptions of the TSP Polytope specifically (see [6], [9], and in particular) are 
not applicable to the developments in this paper. 

Theorem 5 The following hold true: 

1. Ext(Wi) = Wi; 

2. Ext{W L ) = Ext(Wi); 

3. Conv{W L ) = ConviWr). 
Proof. See [TJ or [8]. ■ 

Definition 6 We refer to Conv(Wi) = Conv(Wj) as the "Assignment Polytope." 

As in [3J, in order to reformulate Wj and Wl, we use the framework of the multipartite digraph 
G = (V, A) illustrated in Figure{]\ Each node in this graph is defined by a (city, time-of-travel) 
pair. Arcs are labeled with triplets (a, b, c) G (M, T\{n — 1}, M), and represent the possible choices 
at consecutive times-of-travel. Specifically, the arc linking nodes (i, r) and (j, r + 1) of the graph 
is labeled (i,r,j), and represents the choice to visit cities i and j at travel times r and r + 1, 
respectively. 

Definition 7 

1. The set of all the nodes of Graph G that have a given city index in common is referred to as 
a level of the graph. The i th level (i G M) is denoted L{(G) := {(u, v) G V | u = i}; 

2. The set of all the nodes of Graph G that have a given time-of-travel index in common is 
referred to as a stage of the graph. The j th stage (j £ T) is denoted Sj(G) := {(u,f) £ | 
v = j}; 

3. A path in Graph G that simultaneously spans the levels and the stages of Graph G is referred 
to as a "ciiy and stage spanning (c.a.s.s.) path" of the graph; 

4. The set of all the c.a.s.s. paths of Graph G is denoted f2. That is, 

n := {((ii, 1, i 2 ), (i 2 , 2, i 3 ), (in-2, n - 2, i n _i)) G : i p ^ i q V (p, 9 ) G (T, T\{p})}. 

Remark 8 T/ie association we make in our modeling between Graph G and the TSP is to interpret 
a positive flow into/out of any node of the graph to mean that the corresponding city and time-of- 
travel pair have been assigned to each other. 

Remark 9 

1. There exists a one-to-one correspondence between the perfect matchings of M and T and the 
c.a.s.s. paths of Graph G; 

2. There exists a one-to-one correspondence between the c.a.s.s. paths of Graph G and the 
points of Ext{W L ) = Ext(Wi) = Wi\ 

3. There exists a one-to-one correspondence between the c.a.s.s. paths of Graph G and TSP 
tours. 
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Figure 1: Illustration of Graph G 



2.1 Reformulation of the Assignment Polytope 
Notation 10 : 

1. The set of stages of Graph G from which arcs of the graph originate is denoted R; i.e., 



2. V (i, j,u,v, k,t) G F 6 , V (p, s) G R 2 : 1 < p < s, Zi t \j UVV k s t denotes a non-negative variable 
that represents the amount of flow in Graph G that propagates from arc (i, onto arc 
(k, s, t), via arc (u,p, v); 

3. V (i,j,k,t) e F 4 , V (r, s) G R 2 : 1 < r < s, yi r jkst denotes a non-negative variable that 
represents the total amount of flow in Graph G that propagates from arc (i,r,j) onto arc 
(k,s,t). 

2.1.1 Integer Programming Reformulation 

The constraints of our Integer Programming reformulation of the Assignment Polytope are as 
follows: 



R:=T\{n-l}; 




(3) 



ieM jeM ueM teM 




i, j, k,t,u £ M; p, s£R: l<s<p<n — 2 



(4) 
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^ ] z i,l,jvpukst ^ ] Ju,p+l,ufcst 0) 

i, j, k,t,u £ M; p, s&R: l<p<s — 1, s>3 



Ui,l,jupv ^ ^ ^ ^ z i,l,jupvkst — 0; 

keMteM 

i, j, u, v G M; s€i?:l<p<s 



Ui,l,jkst 

E E = o; 

i,j,k,teM; p,seR:l<p<s 



Ui,l,jkst ^ ^ ^ ^ z i,l,jupvkst ^ ^ ^ ^ z i,l,jkstvpu — 0; 

peR:Kp<sweM p£R:p>sv£M 

i,j,k,t,u£M; s G 



Uupvkst ^ ] ^ ] z i,l,jupvkst — 0; 

v,k,t £ M; p,s G i?, 1 < p < s 



E + E Virjk,r+l,t = 0; i,j<EM; V G i? 

(fc,t)€M 2 (k,t)&M 2 : 

k^j; (k,r+l,t)eA 



E EE^+ E EE^ E E^- 

s&R: s>r k&M teM s&R: s<r keM t£M seR:s>r+l fcgM 

+ E E»+ E E^+ E = 

sei?:s>r+l fceM s£ii:s>r+l fceM s£ii:s>r+2 fceM 

i,j E M; r £ R 

Virjkst^ {0,1}; i,j,k,teM; r,s£R:r<s; 

Zi,i,j U pvkst G {0, 1}; i, j,u, v, k, t e M; p,s £ R: I <p < s 



The propagation of one unit of flow from stage 1 of Graph G is initiated by constraint ([3]). 
Constraints (jl]) and © ensure that all flows initiated at stage 1 propagate onward, to stage n-1 
of the graph, in a connected and balanced manner. Specifically, constraints @ stipulate that the 
total amount of flow from arc that propagates through arc (k,s,t) and enters node (u,p+l) 

is equal to the amount of flow from arc (i, 1, j) that propagates through arc (k, s, t) and leaves node 
1); Constraints © stipulate that the total amount of flow from arc (i, that enters node 
(u,p+l) to propagate on to arc (k,s,t) is equal to the amount of flow from arc that leaves 

node (u,p + 1) to propagate on to arc (k, s, t). Constraints © and ([7]) ensure that the propagation 
of the flow from a given arc at stage 1 of Graph G onto a given arc at another given stage of the 
graph is consistently accounted across all the other stages of the graph. Constraints Q stipulate 
that the total amount of flow that propagates from arc (u,p,v) onto arc (k,s,t) is equal to the 
total of the flows from arcs at stage 1 that propagate onto arc (k,s,t) via arc (u,p,v). Constraints 
([8]) require that the total flow on any given arc of Graph G must propagate on to every level of the 
graph, or be part of a flow propagation that spans the levels of the graph. Constraints (jlOp ensure 
that the initial flow propagation from any given arc occurs in an "unbroken" fashion. Finally, 
constraints (jlip stipulate (in light of the other constraints) that no part of the flow from arc (i, r, j) 
of Graph G can propagate back onto neither level i nor level j of the graph. 

Remark 11 

1. The number of variables in the system (|3j)- (|lip is 0(n 8 ); 

2. The number of constraints in the system (j3j)- (|lip is 0(n ). 

Definition 12 

1. We refer to the set of points in the space of the y- and z- variables that satisfy the system 
(|3|)- (fT3"j) as the "IP Polytope," and denote it by Qj; i.e., Qj := {(y,z) G 7Z> : (y,z) satisfies 
([3j)- (fT3]) }. where £ is the number of variables in the system (f3|)- (fT3|) . 

2. We refer to the linear programming relaxation of Qi as the "LP Polytope," and denote it by 
Ql', i.e., Ql '■= {{y,z) G 7& : (y,z) satisfies (|3]l- (fTT]l . and < (y, z) < 1}, where £ is the 
number of variables in the system ([3|)- (|lip . 

Theorem 13 (y, z) E Qj <^=^> 3 exactly one set of city indices, {i r G M, r = 1, . . . , n — 1}, such 
that: 

1 for r, s G R : r < s, and 

(a,b,c,d)= (i r ,i r+ i,i s ,i s+ i) . , 

Varbcsd = \ (J-4) 

otherwise 



and 



1 for r, s G R : 1 < r < s, and 

(a,b,c,d,e,f) = (ii,i 2 ,i r ,ir+i,is,i s +i) n ^ 

z a,l,bcrdesf — \ \ L0 ) 

otherwise 
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Proof, o) =>: Let (y,z) E Qj. Then, given CE]), and (fT3l) : 
Constraint ([3]) => 3 a unique set of city indices, {i r E M, r = 1, . . . , 4}, such that: 

^11,1,1212.2,1313,3,14 = 1 (16) 
Condition ()15|) follows directly from the combination of (|16p . ((4|), and ([5]). 

Condition (|14p follows from the combination of Condition (|15p with constraints ©-([T]) and ©. 
6) -^=: Trivial. ■ 

Theorem 14 There exists a one-to-one correspondence between the points of Qi and the points of 
Wi. 

Proof. Combining Theorem [13] with constraints (jlip . (|12p . and (|13p . we must have: 

i) ([HI) and (fT5j) are satisfied for (y,z), and 



(y, z) E Qj 3(ii,i 2 , %-i) € M 



n— 1 . 



(17) 

k n) i p / iq V(p, <?) G (T, T\{p}) 



Definition 04 and Condition (I17p ==>There exists a one-to-one correspondence between the points 
in Qj and c.a.s.s. paths of Graph G. The theorem follows directly from the combination of this 
and Remark [9j 2. ■ 



Corollary 15 There exists a one-to-one correspondence between points in Qj and TSP tours. 
Definition 16 

1. We refer to the point w E Wj correponding to (y,z) E Qi as the "assignment corresponding 
to (y, z)," and denote it by the ordered set M(y, z) := (£1, «2, ■ ■ • , i n -i), where i q E M is the 
index of the city visited at time-of-travel q, according to w; 

2. We denote the TSP tour corresponding (y, z) E Qi as T(y,z); i.e., for (y, z) E Qi, T(y,z) 
denotes the travel: 1 — ► i\ — ► 12 — ► • • • — > i n -i — ► 1, where i q E M(y, z) Vg E T. 



2.1.2 Linear Programming Reformulation 

Lemma 17 (Flow conservation lemma 1) Let (y,z) E Q^. The following holds true: V(a, 6) G 
M 2 , V(p, q, r, s) G i? 4 : 1 < p < q; 1 < p < r < s, 

EEEE 

z a,l,b,i p ,p,jp,i q ,q,j q 
i p eM j p <EM i g £M j v £M i r £M j r £M i s GM j s £M 
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Proof. V(a, b) G M 2 , V(p, (j,r,s) 6fi 4 :l<p<g;l<p<r<s, 
E E E E ^,i,fe,v,P,j P ,i 9 ,g,i 9 = E E (Using ®) 

ip£M j p eM i q €M j q €M i p eM j p EM 

= Y Y Y Y Z a,l,b,i P ,p,j P ,ir,r,j r (Using ©) 
i p &M j p £M i r eM j r £M 

= Y Y Y Y z a,i,b,i P ,p,j P ,ir,r,j r (Re-arranging) 

i r <=M j r <=M i p £M j p €M 

= Y Y yaXb,i r ,r,j r (Using ©) 

ir&M j r &M 

= Y Y Y Y z a,l,b,ir,r,jrMj s (Using ©) 
i r <=M j r <=M i 3 £M j a €M 

m 

Lemma 18 (Flow propagation lemma 1) Let (y,z) G Ql. The following holds true: 

V(ii,*2>*3>*4) £ M , 2/11,1,22,13,3,44 > -x=>- 1,12,12,2,13,13, 3,i4 > 0. (18) 
Proof. Using constraints © and (|1U|) . constraints (JT]) for p = 2 and s = 3 can be written as: 

J/ii,i,tai3,3,t4 - ^1,1,1212,2,1313,3,14 = V(«i,i 2 >*3,i4) S M 4 (19) 
The lemma follows directly from (|19p . ■ 

Lemma 19 (Flow propagation lemma 2) Lei (y, z) G Ql. Then, we must have that: 
Vr G R : r > 4, V(n, i 2 , 13, *4, *r, V+l) G 

{0 z n,i,«2,«2, 2, 13, 43,3,14 > 0; 
(20) 
^i) ^j'i, 1,12, i2, 2, i3,i r ,r,i r -\-i ^ 0' 

Proof. 

a) Condition i. Using constraints ©, 

Zi u i,i 2 ,i 3 ,3M,i r ,r,i r+1 > =>■ yii, 1,22,23,3,44 > Vr G i? : r > 4, V(ii, i 2 , i 3 , M, V, V+l) G M% ( 21 ) 
From Lemma [THl 

2/21,1,22,23,3,24 > 2:21,1,22,22,2,23,23,3,24 > V(h,i2,h,u) e M 4 (22) 
Condition i) follows directly from ([22|) . 

6) Condition ii. Using ([U]), (J7|), ([!]), and (fTUj) . constraints (0) for p = 2 and it = is, can be written 
as: 

z h,i,i2,i2,2,i 3 ,i s ,s,is+i - Y z h,i,i2,i 3 ,3,v,i s ,s,i 3+1 = V(«i, i 2 > i 3 , i s , i s+1 ) G M 5 ; Vs G -R : s > 4 (23) 
Hence, in particular, we must have: 

^2l,l,22,22,2,2 3 ,2 r ,r,2 r + l -^2l,l,22,2 3 ,3,24,jr,r,2V + l > V(«l , 2 2 , «3 ) *4 J V, V+l) G ^5 V f G -R i T > 4 (24) 

Condition ii follows directly from (|24p . ■ 
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Lemma 20 (Flow propagation lemma 3) The following holds true for all (y,z) € Ql'- 
Vr € R : 2 < r < n - 4, V(ii, i 2 , V, V+i, V+2, V+3) G M 6 , 

(25) 

M) z ii,l,i2,ir,r,i r +l,ir+i,r+l,i r +2 > 0- 

Proof, a) Condition i. Using ([6]), ([7]), , and (fl~0|) . constraints © for p = r, s = r + 2, and 

u = Vfij can be written as: 

^ ] <2ji,l,i2,t>,r,i r+ i,j r+ 2,r+2,j r+ 3 — ^"i,l,ta,ir+li r+l,i r +2,ir+2,r - +2,i r +3 = ^ 

Vr€i?:2<r<n-4, V(ii, » 2 , Vfl, V+2, V+3) G (26) 
Hence, in particular, we must have: 

z i 1 ,l,i2,ir,r,i r+ i,i r+ 2,r+2,i r+3 ~ z i 1 ,l,i 2 ,ir+i,r+l,ir+2,ir+2,r+2,i r+3 < 

Vr€i?:2<r<n — 4, V(«i,i2> i r , Vfl? V+2, V+3) G M 6 (27) 
Condition i) follows directly from (|27p . 

6) Condition ii. Using constraints ©, ((7J), ©, and (fTU|) . constraints (|3J) for p = r + 1, s = r , and 
u = v+2, can be written as: 

z h,l,i2,ir,r,i r +i,ir+i,r+l,ir+2 ~ ^ ] z ii,l,i2,ir,r,i r +i,i r +2,r+2,v = 

Vr € R:2<r < n - 4, V(ti, » 2) V+l, V+2, V+3) e M 5 (28) 
Hence, in particular, we must have: 

2ii,l,i2,ir,i",i r +i,i r +i,r+l,i r +2 — z ii,l,i2,ir,r,i r +i,i r +2,r+2,i r +3 — 

Vr€i?:2<r<n-4, V(»i, i 2 , V, V+l, V+2, V+3) G M 6 (29) 
Condition ii) follows directly from (f29|) . ■ 

Notation 21 For (y, z) € Ql-' 

1. The sub-graph of G induced by the positive components of (y, z) is denoted as: 

H(y,z): = (P(y,z),E(y,z)), (30) 

where: 
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P(y,z) ■={ 



> U 



(i,r) eV : Y VaXbirj + 

a&M b&M jeM:(i,r,j)eA 

+ Y Y Y Va,l,bj,r-l,i > 

a&M beM jeM:(j,r-l,i)eA 



(31) 



E(y,z):=Ui,l,j)eA:J2yi,ijj,2,t > olu 
I teM J 



i,r,j) eA:^2Y^ Y V^Mrj > 

agM bGM j£M:(i,r,j)eA 

2. The set of arcs of if(y, z) originating at stage r of H(y, z) is denoted T r (y, z); 



(32) 



3. The number of arcs originating at stage r of Graph H(y,z) is denoted j r (y,z) = \T r (y,z) \ . 
For simplicity r y r (y,z) will be henceforth written as 7 r (unless that causes ambiguity); 

4. The index set associated with F r (y,z) is denoted A r (y,z) := {1, 2, . . . , j r }. For simplicity 
A r (y, z) will be henceforth written as A r ; 

5. The v ih arc in T r (y,z) is denoted as a Tjl/ (y,z). For simplicity a r ^(y,z) will be henceforth 
written as a r y, 

6. The tail of a r , v is labeled i r ^ v (y,z); the head of a r ^(y,z) is labeled j r ,v{y,z). For simplicity, 
ir,v(y, z) will be henceforth written as i rji ,, and j r ,u(y, z), as j r ,^; 

7. Where that causes no confusion (and where that is convenient), for r £ iJ with r > 2, and 
(a, p) G (Ai,A r ) "yn iQ ,ij liQ ,i riP ,r,j r , p " will be henceforth written as "y(i,a)(r,p)-" Similarly, for 
(r,a) G i? 2 with 1 < r < s and (a,p,a) G (Ai,A r ,A s ), %, ai ij 1 , a ,i r ,,,rdV, P ,».,a,«d.,/ wil1 be 
henceforth written as a Z(i iQ )(r, / o)(s,o-)-" 

Definition 22 ("Paths in (y,z)") Let(y,z) G Ql. V (r, s) G R 2 : s > max{3, r+1}, V(i/i, i/ r , i/ s ) G 
(Ai,A r ,A s ), a se£ of arcs ofH(y,z), 

{(ar,v r ,...,a StVt ) G (£(y,z)) s - r+1 : z (lil/l)(Pil/p)((?)i , g) > V(p,g) G i? 2 : 
max{2,r} <p<g<s-l; i PjI/p = j p _i, I/p _ 1 Vp G (-R n [r + l,s])} 

is referred to as a "path in (y,z) from (r,v r ) to (s,v s )." 

Notation 23 Let (y, z) G Q L - V (r, s) G R 2 : s > max{3, r + 1}, V (p, a) G (A r , A s ), 

1. The set of all paths in (y,z) from (r, p) to (s,cr) is denoted U^ rp ^ scr ^(y, z); 

2. The index set associated with Ut rp \t sa \(y, z) is denoted ^( ri p)( S)(T )(y, 2) := 
{1,2,..., <P(r,p){ s ,c) (y, «)}, wh ere <P( r , P )(s,a) (y, «) := | ^(r, P )(s,a) (y, *) I ; 
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3. The k th element of U( r , P )(s,a)(y, z) G $( r ,p)(s,a){y, z)) is denoted C(r tP ^ s ^ }k (y, z). 

Theorem 24 (Distinct paths in (y,z)) Let (y,z) G Ql- Then, V (r, s) G R 2 : s > max{3,r + 
1}, V(ai,a 2 ) G A 2 , V(/3\,/3 2 ) G A 2 : (£/ (r , ai)Ml) (2/, z) + and U (rta2)(t>M (y, z) ^ 0), V k 

{5 G R : r < 5 < s; 



(7i,72) ^ (A 9 ,A 9 \{7 1 })} 9 |a 9i7i G £ (r . iai)i(Si/3l)ifc (y,2;) and a 9i7a G £( r ,a 2 ),( s ,/3 2 ),t(y, *)J • 

Proof. The theorem follows directly from the combination of constraints ([6]) , (JT]) , ([9]) , and (|10p , 
and Definition [22j ■ 

Theorem 25 (Path structure theorem 1) Let (y,z) G Q^. The following holds true: 
V(r,s) G i? 2 fflii s > max{3,r + 1}, V(p,cr) G (A r , A s ), y( r , p )( s , a ) > #(r,p)(s,<r)(y> z) / 0. 

Proof. First, note that it follows directly from the combination of Lemmas 1181 [191 and [201 that 
the theorem holds true for all (r, s) G -R 2 with s G {r + 1, r + 2}, and all (u r , v s ) G (A r , A s ). 

a) =x Assume there exists an integer uj > 2 such that the theorem holds true for all (p, t) 
G R 2 with t = p + uj, and all (u p , vt) G (A p , At). We will show that the theorem must hold for all 
(p, u) G R 2 with u = t + 1 = p + o> + 1, and all (v p , v u ) G (A r , A u ). 

Let (p, G -R 2 with u = p + cj + 1, and (u p , u u ) G (A p , A u ) be such that: 

y( P ,u p )(u,u u ) > 0. (33) 



Define: 



B( P ,v p )(u,u u )(y,z) := {a G Ai : Z(i^){ P , Vp )(u,v u ) > 0}. 



(34) 



Then, Q and (pi) 



wi t h 



< 




(35) 



y( P ,u p )(u,u u ) 



z (l,a)(p,v p )(u,v u ) 



Condition ([35j), and constraints ([5]) and fp~Q]> 



*) ip+1,/3 = j P ,v P V /? G C a) ( P) ^)( u ^)(?/,^) 



)(y,z); and 



(36) 



m J 2(l,a)(p,K p )(u,l/„) 



< 




^(l,a)(p+l,/3)(«,i/„)- 
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By assumption (since u = (p + 1) + u), condition (f36U i) => 

)(v,z) (37) 
Also, it follows from the combination of condition (|33|) . constraints ©, and constraints ([5]), that: 

V/?eC«,( w )( U)1 ,j(t(,4 3 {Y(p, Vp )(p+i,p)(u,vu)(y, z ) £®(p+iMwu)(y> z> )} 3 : 

{( z )(l, a )(p,u p )(q,u q ,,) > Vi G T(p )I/p )(p +1)j g)( W)l/tt )(2/,2), 

Vg G (i? n [p + 1, it]), and V i/ 3jt G A 9 : Op,^^ G £(p+i,/3)( U)1 , u ), t (z/, z)} (38) 

Hence, V /5 G C a) (p )I/}) )( U)1/tt )(y, z) and V t G T(p,u p )(p+i,p)(u,u u )(y, z ), 

L := (£ {p+m{UyUu))i (y,z) U {%,,„„}) 

is a paf/i in (y ; ,z) from (p, v p ) to (u,u u ). Hence, we have that ^( Pji / p )( UjI / u )(y, z) ^ 0. 
b) <=: Follows directly from Definition (|22p and constraints ©. ■ 

Corollary 26 Fei (y, z) G Ql- The following hold true: 

i) V s G #\{1} , V(a,cr) G (Ai,A s ), y(l, a )( S)0 -) > ^(i,a)(s,<r)(y^) 7^ 0- 

iij V (r, s) G (i?\{l}) 2 with s > max{3,r + 1}, V(a,p,cr) G (Ai, A r , A s ), 

{ii.l) U^ a ^ s ^(y,z) ^ 0, and 
u.2) 3 « g *(i, a )( S)(T )(y,^) 3 ar !() e£(i,a),( ¥ ),j!/,4 

Definition 27 ("TSP tour in (y,z)") Let (y,z) G Ql. V(^i,i/ n _ 2 ) G (Ai,A n _ 2 ), a path in 
(y, z) from (1, u{) to (n — 2, ^-2) is referred to as a "TSP tour in (y, z) (from (1, v\) to (n, i^ n - 2 )). 77 

Notation 28 Let (y,z) G Ql- For all (a, (3) G (Ai,A n _ 2 ), 

1. The set of all paths in (y, z) from (1, a) to (n — 2, 0) is denoted as n Q ^(y, z); 

2. The index set associated with IT a/ g(y, 2) is denoted ^f a p(y,z) := {1, 2, ... , 7T a p(y, z)}, where 
■K a p{y,z) ■= \U a p(y,z)\; 

3. The element of H a p(y, z) is denoted V a /3k(y, z). 
Remark 29 Let (y,z) G Q L - V (a,/?) G (Ai,A n _ 2 ), 

i- n Q/3 (y,z) = c/( lia )( n „2,/3)(y,2); 

2. ^ a p{y,z) = *(i,a),( w -2^)(y,^); 

3. n a/ 3(y,z) = y(i )a )( n _2„g)(y,z); 

4. We assume (w.l.o.g.) that: V af 3k{y,z) = £(i, a ),( n -2,i3),k(y> z ) V ^ G Vap(v, z )- 
Theorem 30 (Equivalences for TSP tours in (y,z)) For (y, z) G Ql : 
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i) every TSP tour in (y, z) corresponds to exactly one perfect matching of M and T; 

ii) every TSP tour in (y, z) corresponds to exactly one TSP tour. 

Proof. Definition 1221 constraints (jlip . and Remark [HJ and Definitions [71 3-d 4 imply that every 
TSP tour in (y, z) corresponds to exactly one c.a.s.s. path of Graph G. The theorem follows directly 
from the combination of this and Remark [9j ■ 

Theorem 31 ("Convex independence" of TSP tours in (y,z)) Let (y,z) G Ql. A given TSP 
tour in (y,z) cannot be represented as a convex combination of other TSP tours in (y,z). 

Proof. Theorem 1301 implies that every TSP tour in (y,z) coresponds to an extreme point of 
Conv(WL) = Conv(Wj). The theorem follows directly from this. ■ 

Theorem 32 (Path structure theorem 2) Let (y,z) G Ql- The following holds true: V r £ R, 
VpG A r; 3 {(a,/3) G (Ai,A n _ 2 ); i>eV a p(y,z)} 3 {a r>p G V a p L (y, z)} . 

Proof. 

Case 1: r = 1. From (El and <Mh: 



a r , P G E(y, z) => 3 a G A 2 3 y( r ,p){2,a) > 0- (39) 
Condition (|39p and constraints ([6]) ==>■ 

3 7 S A n _ 2 3 2(r,p)(2, Q )(n-2, 7 ) > 0. (40) 

Condition (|40p and constraints ([6]) ==> 

3 7S A n _ 2 3 y( r ,p)(n-2,j) > 0. (41) 
The theorem follows from the combination of (|4ip with Theorem 1251 
Case 2: r = n-2. From (JEE]) and ([32]): 

a r , P G £(y, z)=^3a G Ai 3 y(i, a )(r,p) > °- ( 42 ) 
The theorem follows from the combination of (|42l) with Theorem 1251 
Case 3: 1 < r < n-2. From ([3TJ) and (|32|): 

Or,p G E(y, z) 3 a G Ai 3 y(i, a )( r , p ) > 0. (43) 
Condition (f4l3j) and constraints ([6]) 

37 G A n _ 2 3 Z(l, a )(r,p)(n-2, 7 ) > 0. (44) 
The theorem follows from the combination of (|44h with Corollary \2Qli. ■ 
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Corollary 33 Let (y,z) G Ql- The following hold true: 
i) Va G Ai, Vr G R\{1}, Vp G A r , 

2/(i,a)(r,p) > "<=>• 3 (f3 G A n _ 2 ; *■ G ^ a p(y,z)) 9 a r , p G V a f3L(y,z); (45) 
tt) V(r, s) G i? 2 : 1 < r < s, V(a, p, a) G (Ai, A r , A s ), 

z (i,a)(r,p)( s ,a) > ^> 3 (/? G A n _ 2 ; i G ^ Q/3 (y, z)) 3 (a r , p , a S)CT ) G T^G/, z )- ( 46 ) 
Lemma 34 (Flow conservation lemma 2) Lei (y,z) G Ql- The following hold true: 

*) y(l,a)(r,i/ P ) = 2 S 

fpGAp /3gA„_ 2 te*a/3(j/,«): 

Va G Ai, V(p, r) £ E 2 : 1 < p < r, V z/ r G A r 

") V(l,a)(r,iv) = 2 (l,a)(r,^)( g ,^) 
vq&hq /3eA„_2 te* a/3 (y,2;): 

Va G Ai, V(r, g) G i? 2 : 1 < r < g, V v r G A r 

Proof. The lemma follows directly from the combination of constraints ([6]) and ([7]), and 
Theorems [Ml EH and E21 and Corollary E3 ■ 

Definition 35 ("Weights" of TSP tours in (y,z)) Let (y,z) G Q L - For (a, (3) G (Ai,A n _ 2 ) 
such that yn ia )( n -2,p) > 0, and k G ^?i,n-2(y, z), we refer to the quantity 

uj a p k {y,z) := mm {z(i, a )(p,u p )( q ,v q )} (47) 

(p,q)£R z ; {v p ,uq) £ (Ap,A q ): 
l<P<9i (Op,^: a,,i/ g ) £ Plp k {y,z) 

as the "weight" of (TSP tour in (y,z)) Va/3k(y,z)- 

Remark 36 It follows directly from Definitions [2B and\35\ that for (y,z) G Ql, ^ a /3i(y,z) > V 
{a, /3) G (Ai,A n _ 2 ) : Vapfaz) ^ 0, V t G P a p L {y,z). 

Theorem 37 (Path structure theorem 3) Xei" (y, z) G Ql. JTie following hold true: 
i) V a G Ai, V r G i2\{l}, V p G A r , 

y(i,a)(r,p) = X] w a/3*(l/,*) 



"r, P e V aBL (y,z) 



ii) V(r, s) G i? : 1 < r < s, V(a, p, a) G (Ai, A r , A s ), 

S(l,a)(r,/0(s,<r) = U af3lj (y,z) 
/3eA n _ 2 te* Q/ 3(y,2): 
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Hi) V(r, s) G R 2 : 1 < r < s, V(p, cr) G (A r , A s ), 

y(r,p)(s,a) = Y Y Y U a p L (y, Z ) 

(ar, p , a a , a ) £ V 2 afiL {y,z) 

Proof, a) Condition i. First, note that from the combination of constraints ([3]), (JJJ), ([5]), and 
([10]); Remark [36] and Theorems (f3Tj) . and ([32]) . we must have: 



E E E «(l,a)(2,/3)(3,«) = Y Y w a fit (y,«) = l- 

aGAi /3eA 2 : <5eA 3 : a£Ai p£A n _ 2 te*a e (l/,*) 

h,0=h,a h,S=32,p 



(48) 



a.l) From Definition 



z (i,")(2,/3)(3,5) = X] Y u aQi (y,z) VaGAi 

/3gA 2 : <5eA 3 : gGA n _2 tG* a£J (j/,2) 



(49) 



Lemma [T7] and relations (|49p ==>■ 

X *(l,«)(r,p)(s,o-) = Y u agL(y,z) Va G A x , V(r, s) G i? 2 : 1 < r < s (50) 

peA r aeA s £>gA„_ 2 t,em ae {y,z) 

Using constraints ©, relations ([50j) =>■: 

S ^.oJM = S Y ^(y^) VaG Ai,Vr Gi?\{l,n-2} (51) 

peA r £>£A n _ 2 u&S ae {y,z) 

Using Theorem 1241 (|5ip can be written as: 



^ V{\,a){r, P ) = Y Y u agL {y,z) Va G Ai, Vr G -R\{1, n — 2} 

peA r p£A r £>£A n _ 2 te* ae (2/,«): 

a r ,p £ VnQi.{y,z) 



(52) 



Re- arranging (|52p gives: 



E 

p£A r 



\ 



y(l,o)(r,p)- X] zE w Qet (y,z) 



V 



a r , P £ Va B L{y,z) 



Va G Ai,Vr G 22\{l,n - 2} (53) 



/ 



a. 2) Combining Lemma [341 n with Definition 1351 we have that: 

l/(l,a)(iy>) = Y Y Y Z (l,*)(r,p)(q,v q ) 

v q eA q £i£A„_2 i£^ ae (y,z): 

>- E E E 

^ 9 £A«j £>£A„_2 t£*a e (j/,«): 

(ar, P , a,,,,,) £ P^Cs/,*) 



> XI zE 

£>£A n _ 2 tG* ae (j/,2): 

Or,p £ Vagi,{y,z) 

V r G n - 2}, Vq e R : q > r, V(a, p) G (Ai, A r 



(54) 
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Relations (1531) and (15 



y(i,a)(r, P )= ^2 ^2 u aBl ,(y,z) Vr g R\{1, n - 2}, V(a, p) g (Ai, A r ) 

£>eA n _ 2 ie* ae (y,2;): 
a. 3) Using constraints ([7]), relations ([50]) =>: 



y(i,a)( s ,o-) = X] X] ^agiiy^) Va € Ai , Vs € i? : s > 2 
o-eA s £>eA„_2 ie* ae (y,z) 

Using Theorem [Ml ([56]) 

= X] X! u agL (y,z) Va € Ai,Vs G i? : s > 2 

°-eA s o-eA s e€A n _2 i6* ae (?/,z): 

as,tre'Pa()t(j/,2) 



Re-arranging ([57]) gives: 

/ 



E 

ctSA s 



\ 



a s ,cr6'P a e t (j'!' z ) 



Va G Ai,Vs € i? : s > 2 



/ 



o.4) Combining Lemma IH^li with Definition 1351 we have that: 

= E E E 

^ P eAp £ieA n _2 te*a e (l/,«): 

(ip,„ p G V al 3 L (y,z) 

aEE E 

> X] £ u agL (y,z) 

£>eA„_ 2 j.e* afi (2/,z): 

V(p, s) G i? 2 : 1 < p < s; s > 2, V(a, cr) G (Ai, A s ) 
o.5) Relations ([58]) and (JSSJ) => 



= X] u aeb (y,z) Va G Ai,Vs G -R : s > 2 

eeA„_ 2 (.e* ae (y,z): 

a s ,<r€.Va.Q<.{y,z) 

a. 6) Condition i of the theorem follows from the combination of (|55p and (|60p . 

6) Condition ii. 
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6.1) Using Theorem 1241 and Corollary 1331m. (|50p can be re-written as: 

E E Z E "afr{y,z) = E z (l,a)(r,p)(s,v) 

peA r creA s /3gA„_ 2 te* Q/ 3(j/,z): peA r o-eA s 

V(r, s)eE 2 :l<r<s, Va G A x 
Re-arranging ([6"Tj) gives: 
/ 



(61) 



EE 

p£A r <tGA s 



^(l,a)(r,p)(s,(7) - E 



V 



^2 u a po(y,z. 



V(r, s) G R 2 : 1 < r < s, Va G Ai 

6.2) Prom Lemma [34"1m. we have: 

y(l,a){s,a) = z (l,a)(r,p){s,a) + E E E z (l,a)(r,u r )(s,a) 

u r eA r -. /3eA„_ 2 t.e^ al3 (y,z)-. 

V(r, s) G i? 2 : 1 < r < s, V(a,p,cr) G (A x , A r , A s ). 

6.3) From Condition i, we have: 



(62) 



V(l,a)M = E E 



00 



a/3. 



EE E 



Vr+P (a r ,, r ,a a ^)eP^ L (y,z) 



V(r,s) £ fi 2 : 1 < r < s, V(a,p,a) G (Ai,A r ,A s ) 

6.4) Definition [35] 

EE E 2 E E 

v r &A r : /3GA n _2 te* a/3 (j/,2;): 

(ar,^,a., ff )eJ^ /3t (y,*) 

V(r,s) G i? 2 : 1< r < s, V(a,p,a) G (Ai,A r ,A s ) 

6.5) Relations f[H5|t - fjB5|> 

Z(l,a)(r,p)(s,a) < E E U af3L (y, Z ) 

/3eA„_ 2 te* a/ s(l/,«): 

V(r,s) G -R 2 : 1 < r < s, V(a,p,a) G (Ai,A r ,A s ) 

6.6) Combining (|62p and (|66p . we must have: 

z (l,a)(r,p)(s,u s ) = E E U al3i (y,z) 

/3eA n _ 2 ief a/ 3b,z): 

(or, P ,o s ,o-)eP^ i g i ,(y,z) 

V(r,s) 6 U 2 : 1 < r < s, V(a,p,cr) G (Ai,A r ,A s 



E 



(a r ,„ r ,a a , a )£Plp L {y,z) 



(63) 

(64) 

w a /3 t (y,z) 
(65) 



(66) 



(67) 
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c) Condition Hi. From the combination of constraints (O) and Condition ii, we have: 



y(r,p)(s,a) = Z(l,a)(r,p)(s,a) = Yl E E U a p L (y,z) 

(or, p , a s , v )eP^ 0l (y,z) 

V(r, s) G R 2 : 1< r < s, V(p, a) G (A r , A s ). 



Corollary 38 (y, z) G Ql (j/)^) corresponds to a convex combination of perfect matchings of 
M and T with coefficients equal to the weights of the corresponding TSP tours in (y,z). 

Theorem 39 The following holds true: Conv{Qi) = Conv(Qj). 

Proof. The theorem follows directly from the combination of Theorem [30| Theorem [31~j and 
Corollary [ 



Corollary 40 The following mappings are bijective: 



1. Bx : Conv{Q L ) 

2. B 2 : Conv(Q L ) 

3. B 3 : Ext(Q L ) h 

4. £ 4 : Ext(Q L ) h 

5. £ 5 : Ext{Q L ) h 

6. £ 6 : Ext{Q L ) h 



-» Com;(Wj); 
-» Com;(iy L ); 
Ea^Com^W})); 

O; 
A. 



2.2 Reformulation of the Travel Costs 

Definition 41 (Re-defined costs ) Let (y,z) G Ql- V(i,j,u,v,k,t) G M 6 , V(r, s) £ R 2 : 1 < 
r < s, the "cost" associated with z u \^i r jk s t is defined as: 

%l, u ~\~ %u,v ~\~ %v,j ~\~ %j,t if r = 2; s = 3; i = vj k = j 

Xjj if r = 2; 4 < s < n — 3; £ = u 

Xj.t + ^t.i if r = 2; s = n — 2; z = w 

Otherwise 
Theorem 42 Lei: 



Cul,virjkst 



l,virjkstZu,l,virjkst' 



virjk 



u£M v£M i€M reR: jeM keM seR: t£M 
r>l s>r 



Then, for (y,z) G Ext(Qi,), accurately accounts the cost of the TSP corresponding to (y, 

z). 
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Proof. From Theorem 1391 

(y, z) G Ext(Q L ) ^ (y, z) E Qi (70) 
Now, using Theorems I13| it can be verified directly that for (y,z) G Qj, 

n-2 

x i r ,i r +i "I" ^-in—ijlj where i r G AA(y,z) V r G T. ■ 

r=l 

2.3 Overall Linear Program 

Our overall linear programming model is as follows: 

Problem 43 (Problem LP) 

min {ft(y,z) : (y,z) G Q L } 

Theorem 44 The following statements are true of basic feasible solutions (BFS) of Problem LP 
and TSP tours: 

1. Every BFS of Problem LP corresponds to a TSP tour; 

2. Every TSP tour corresponds to a BFS of Problem LP; 

3. The mapping of BFS's of Problem LP onto TSP tours is surjective. 

Proof. Statements (1441 1) and (1441 2) follow directly from the combination of Theorem [39| 
Corollary [3016, and the correspondence between BFS's of LP models and extreme points of their 
associated polyhedra (see [1] pp. 92-101]). Statement (f4"4l 3) follows from the primal degeneracy of 
Problem LP (see [gj p. 32]). ■ 

Corollary 45 Problem LP solves the TSP. 

3 Numerical Implementation 

We implemented a streamlined version of Problem LP in which constraints ([9])- (|ll|) were handled 
implicitly, and the upper bounds on the variables were omitted. We solved 3 problems with n = 8 
cities, and 3 problems with n = 9 cities. One of the 8-city problems and one of the 9-city problems 
had, each, travel costs all equal to zero (i.e., we had ty = V(«,j) G (N, N\{i} in these problems). 
Each of the remaining four of our test problems had travel costs that were randomly generated 
between 1 and 300. Also, for each problem size, one of the randomly-generated problems had 
symmetric travel costs, while the other had asymmetric costs. 

We used the simplex implementation of the "Clp" routines of the COIN-OR open source library 
([2]) to solve the primal and dual forms of each of the test problems, respectively. The computational 
results are summarized in Table 1. In general, the dual forms performed better than the primal 
forms. The computational times across the two forms range from 1.02 seconds to 34.50 seconds 
(Sony VAIO notebook computer with 1.83 GHz Intel Core 2 Duo processor) for the 8-city problems, 
and from 33.84 seconds to 4,818.44 seconds for the 9-city problems. Comparable times are on the 
order of several hours (more than 6, in general) for 8-city problems for the model in [3] . 
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J. LUUlvlll 


l\Fi I'm nr»T" 


Tvr>p nf 


Primal Form 


Dual Form 


TMlirnHpr 


of pi tips 


cost a 


^ of iterations^ 


CPU seconds c 


^ of iterations^ 


CPU seconds 


1 


8 


A 


15,354 


11.64 


1,274 


1.02 


2 


8 


S 


11,227 


9.23 


443 


0.56 


3 


8 


N 


10,265 


2.84 


10,610 


34.50 


4 


9 


A 


149,257 


2,653.92 


3,145 


6.77 


5 


9 


S 


217,914 


4,818.44 


70,598 


384.25 


6 


9 


N 


32,831 


33.84 


69,498 


2,149.55 



a: "A" = asymmetric; "S" = symmetric; "N" = all costs are equal to zero 
b: Number of iterations 

c: Sony VAIO notebook computer (1.83 GHz Intel Core 2 Duo processor) 



Table 1: Summary of the computational results 

4 Conclusions 

We have developed a new linear programming model of the TSP. The formulation is an order of 
size smaller than that of the existing model, and results in computational times that are several 
orders of magintude smaller. Hence, we believe the proposed formulation represents a move in the 
right direction with respect to eventually being able to solve more realistically-sized TSP's. 
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